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ABSTRACT 

We have used the Gemini Multi-object Spectrograph (CMOS) on the Gemini 
South 8m telescope in exceptional conditions (0.6" FWHM seeing) to observe the 
outer stellar disk of the Sculptor group galaxy NGG 300 at two locations. At 
our point source detection threshold of r' = 27.0 (3ct) mag, we trace the stellar 
disk out to a radius of 24', or 2.2 R 25 where R 25 is the 25 mag arcsec”^ isophotal 
radius. This corresponds to about 10 scale lengths in this low-luminosity spiral 
(M^ = —18.6), or about 14.4 kpc at a cepheid distance of 2.0 ± 0.07 Mpc. The 
background galaxy counts are derived in the outermost held, and these are within 
10% of the mean survey counts from both Hubble Deep Fields. The luminosity 
prohle is well described by a nucleus plus a simple exponential prohle out to 10 
optical scale lengths. We reach an effective surface brightness of pip = 30.5 mag 
arcsec”^ (2cr) at 55% completeness which doubles the known radial extent of the 


Fig. 1.— Wide held image (30'x20') of NGG 300 from the Digitized Sky Survey (North is 
at the top and East is left). The insets show the summed data from the two GMOS outer 
helds. The GMOS held of view is 5.'5. The helds are located at a radius of 15.47' (Field 1) 
and 21.14' (Field 2). 
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optical disk. These levels are exceedingly faint in the sense that the equivalent 
surface brightness in i? or K is about 32 mag arcsec”^. We hnd no evidence 
for truncation of the stellar disk. Only star counts can be used to reliably trace 
the disk to such faint levels, since surface photometry is ultimately limited by 
nonstellar sources of radiation. In the Appendix, we derive the expected surface 
brightness of one such source: dust scattering of starlight in the outer disk. 

Subject headings: galaxies: individual (NGC 300) — galaxies: photometry — 
galaxies: formation, evolution, dynamics, structure — ISM: dust 


1. Introduction 

The structure of the outer parts of galactic disks is central to our understanding the 
formation and evolution of spiral galaxies. The surface brightness profile of a spiral galaxy 
is a hrst order signature of the formation process (Freeman & Bland-Hawthorn 2002; van 
der Kruit 1987; Ferguson & Clarke 2001). By measuring scale lengths for galaxy samples 
of different redshifts, it is for example possible to directly address the evolution of galactic 
disks (Labbe et al. 2003). 

To a hrst approximation, the radial light distribution of disks is well described by an 
exponential decline (Freeman 1970). If the exponential is due to dynamical evolution (Lin 
& Pringle 1987), it is important to establish out to what radius this process took place. It is 
possible that departures from an exponential will show up at faint levels, rehecting new forms 
of dynamical evolution or different angular momentum of the outermost accreted material. 

Previous studies have shown that the exponential light distribution breaks smoothly at 
a well-dehned ‘break radius’ in some galaxies (Pohlen et al. 2002; de Grijs et al. 2001). In 
the optical, this occurs at galactocentric distances of typically 3-5 exponential scale lengths. 
In principle, this break radius could be used as an intrinsic parameter measuring the size of 
a galactic disk. However, despite the fact that these breaks have been known for more than 
20 years (van der Kruit 1979), we do not have a secure understanding for the origin of the 
break, so their application as a fundamental disk characteristic is therefore not yet possible. 

Beyond the break radius, stars are detected albeit with a steeply declining luminosity 
prohle {e.g. Kregel, van der Kruit & de Grijs 2002). To date, the outer edge of the stellar 
disk (‘truncation radius’) has not been detected in any galaxy, largely because diffuse light 
photometry is unreliable at such faint magnitude limits. We note that in two previous papers, 
optical disks are seen to extend to 10 optical scale lengths. In their study of NGG 5383, 
Barton & Thompson (1997) reach down to almost /xy ~ 30 mag arcsec“^ (1<^)- In a multiband 



3 


study of NGC 4123, Weiner et al. (2001) reach 29, 29, 28 mag arcsec“^ (la) in B, V and I 
respectively. However, we point out that both of these galaxies are strongly barred and have 
peculiar morphologies. 

We choose to focus our hrst study of outer disks on the Sculptor group for several reasons. 
It is the closest galaxy group to the Local Group with a high proportion of disk systems. 
Furthermore, the group has not undergone much dynamical evolution — the galaxies are 
more or less isolated. Their dark haloes are thought to extend far beyond the bright disks. 
In principle, accretion onto the halo could give rise to disks which are much larger than 
normally assumed. 

The layout of the paper is as follows. Our motivation for targetting NGG 300 is discussed 
in §2. We briefly discuss surface photometry vs. star counts in §3 and show that the 
latter is the only way to get below 30 mag arcsec“^ reliably, an issue we return to in the 
Appendix. The observations and reductions under IRAF are discussed in §4, while details 
of the photometric calibration are given in §5. The results are presented in §6 and these are 
discussed in 57. 


2. Object selection 

We chose NGG 300 as a relatively undisturbed member of the closest galaxy group. 
Sculptor. NGG 300 lies at the nearside of a distribution which is elongated along the line 
of sight (Jerjen, Freeman & Binggeli 1998). There have been various distance estimates 
based on cepheid variables (Freedman et al. 2001; Gieren et al. 2004) and the tip of the red 
giant branch (Butler et al. 2004). We adopt a weighted mean of 26.5±0.1 for the distance 
modulus, equivalent to 2.0±0.1 Mpc. We show below that our photometry reaches down to 
an intrinsic stellar magnitude of M^. = 0.5 (3a) which includes some of the giant branch and 
much of the horizontal branch in a moderately metal poor ([Fe/H]«i-1) population. 

The late type spiral galaxy NGG 300 is ideal for our proposed study. It is an almost pure 
disk galaxy (bulge light < 2%) and has a mild inclination [i = 42°). The galaxy lies at high 
galactic latitude (6 ~ —79°), and therefore has low foreground reddening of E(B-V)=0.013 
mag. NGG 300 shares important similarities with M33 in the northern hemisphere (Blair & 
Long 1997). The optical disk scale length has been measured as 2.17 kpc in Bj (Garignan 
1985; scaled to a distance of 2.0 Mpc) and 1.47 kpc in / (Kim et al. 2004). We note that 
these scale lengths would indicate that the colour of the outer disk is likely to be very blue, 
but this has not been measured to date (see §7). Garignan’s (1985; Fig. 8) photometric study 
reveals that the outermost contours are well behaved, with minimal asymmetry. 
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3. Surface photometry vs. star counts 

There are two basic methods for detecting the faint outer reaches of galaxies: deep 
diffuse light imaging {e.g. de Vaucouleurs & Capaccioli 1979), and faint star counts {e.g. 
Pritchet & van den Bergh 1994). 

Surface photometry of face-on galaxies provides radial information which is more easily 
interpreted, but is limited by the brightness of the night sky. Edge-on galaxies are easiest to 
search for optical breaks and truncations because of the enhanced surface brightness from 
the line-of-sight integration. However, with edge-on systems one cannot be sure that the 
surface photometry is not affected by smaller-scale structures in the disk, like faint spiral 
arms for example, which might mimic the effects of truncation (Pohlen et al. 2004). Dust 
extinction is also more severe in edge-on systems although less so at large radius. 

While diffuse light imaging is a direct method, there are many technical difficulties 
at levels below hr ~ 27 mag arcsec”^ {e.g. Pohlen et al. 2002). These include scattered 
light internal to the optics, difficulties of flathelding data (Pohlen et al. 2002), zodiacal light 
(Bernstein et al. 2002), and diffuse Galactic light (Haikala & Mattila 1995), a combination of 
dust-scattered starlight, dust luminescence, and emission from cool gas. Over the sky, there 
are also variations in Galactic cirrus (Schlegel, Finkbeiner & Davis 2000), and warm ionized 
gas due to the Reynolds layer (Reynolds 1992). At fainter levels, in the range fiR ^ 30 — 32 
mag arcsec“^, the outer HI envelope may have (i) a sufficient column of dust to scatter 
detectable levels of starlight from the optical disk (see Appendix), or (ii) sufficient levels of 
ionization due to a local or global radiation held to produce huorescent emission (Maloney 
1993; Maloney & Bland-Hawthorn 2001). 


A more powerful method is to count individual stars in the old stellar population. The 
technique is particularly effective since good seeing allows contamination from background 
galaxies and foreground stars to be properly evaluated. In their study of the M31 halo 
{{m — M)o ~ 24.5), Pritchet & van den Bergh (1994) achieved 50% completeness for point 
sources down to V = 24.5 in 1" seeing. This corresponds to an ‘effective’ surface brightness 

Fig. 2.— Rejection masks for the GMOS helds. The elliptic arcs illustrate how the star 
counts are binned in azimuth. The major axis of the annuli passes through the middle of 
Field 1. For consistency, all radial prohles are presented in terms of the major axis distance 
from the galactic nucleus {i.e. dehned by the major axis radius of the elliptic annulus). The 
region eclipsed by the on-instrument wavefront sensor can be seen in the top right hand 
corner of both helds. The {x, y) axes are in units of GGD pixels. 
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limit of nv = 29 mag arcsec”^ after correcting for foreground contamination and the stellar 
disk contribution from unresolved stars. With their 4m ground-based observations, Pritchet 
& van den Bergh reach to My ~ 0, corresponding to roughly 40% of starlight from the 
old stellar population. More recently, very deep Hubble Space Telescope observations reach 
mi = 31 at 50% completeness in a relatively small 3' held (Brown et al. 2003). 

With stellar photometry, we can trace the disk to much fainter effective surface bright¬ 
nesses by resolving individual stars. This is the motivation for our deep CMOS observations. 
It is however crucial to establish the magnitude completeness limit and the fraction of missing 
light below the detection limit, which we discuss. 


4. Observations and Data Reduction 

Deep r' images of NGC 300 were obtained in two outer locations using Gemini Multi- 
Object Spectrograph (GMOS) on Gemini South during semester 2003B. The two GMOS 
helds lie close to the major axis of NGG 300 and span a major axis distance range of 7 
to 16 kpc. Our science goal was to track the radial prohle of the galaxy as far as possible 
which required the best possible observing conditions. Of the available hlter set, the r' hlter 
(Fukugita et al. 1996) is the most efficient since (i) this is where the performance of GMOS 
peaks (Hook et al. 2004), (ii) the resolved stellar population will probably be dominated by 
red giants, and (iii) the sky in r' is signihcantly darker than in redder bands like F. 


The observations are summarised in Table 1. The field locations are illustrated in 
Fig. 1. The inner field (Field 1) was placed at the optical edge identihed in the Digitized Sky 
Survey (DSS) to provide continuity with earlier photographic work. While the DSS limit is 
about 25.5 mag arcsec“^ in Bj (3(t; Gorwin 1980), the digitized survey can be displayed at 
high contrast to reveal a more extended disk. The disk edge is broadly consistent with the 
photographic limit from photographic amplihcation (DF Malin, personal communication). 


Fig. 3.— A 25" X 25" subset from the centre of Field 1 to demonstrate the quality of the 
daophot photometric analysis. The images show (a) before source subtraction; (b) with 5a 
sources subtracted; (c) with 3a sources subtracted. The psf is mildly elliptic (10 — 15%) and 
0.”60 along the long axis. 
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The outer field (Field 2) was made to overlap slightly with Field 1 in order to check the 
relative photometric calibration. 

The CMOS detector comprises three 2048 x 4608 EEV CCDs arranged in a row where 
the long dimensions are edge butted together (Hook et al. 2004). After 2x2 on-chip binning, 
the pixels are 0.146" in size. The square field of view in a single CMOS image is 5.5' on a side. 
The on-instrument wavefront sensor and small CCD pixels provided a corrected psf of 0.60" 
FWHM for Field 1, and 0.63" FWHM for Field 2. The psf was found to be mildly elliptic at 
about the 10—15% level. A small psf is essential for removing contaminating sources in the 
field. The background levels between exposures varied by no more than 5% generally and 
the photometric conditions were excellent. 

Each of the fields was observed in 9x900 sec dithered exposures. This allowed us to 
remove cosmic ray events and edge effects arising from the three edge-butted CCDs. Initially, 
we used the standard IRAF/Gemini reduction package. The processing pipeline includes bias 
subtraction and fiat fielding (gireduce), mosaicking of individnal CMOS CCDs into a single 
reference frame (gmosaic), followed by combining the 9 dithered exposnres (imcoadd). The 
images were then trimmed to the area of common coverage. The final snmmed images are 
shown as insets in Fig. 1. 

The FITS data files arrived withont field position information in the headers. Onr 
next step was to astrometrically calibrate the snmmed field images. This was made possible 
by using the GAIA package developed by Peter Draper at the University of Dnrham. We 
identified 20 stars in each field with known astrometric positions in the United States Naval 
Observatory (USNO) photographic catalogne (Monet et al. 2003). 

Once the fields were correctly aligned with respect to each other, we formed a mask for 
each field in order to reject nnwanted pixels. Both masks are shown in Fig. 2. We remove 




a (J2000) 

S (J2000) 

Date (2003) 

Exposnre (s 

Field 1 


00''55™ 50" 

-37° 47' 20" 

Dec 20-22 

9 X 900 

Field 2 


00''56™15" 

-37° 46' 49" 

Dec 23-26 

9 X 900 

Landolt 

95 

OO'' 53™ 01" 

+00° 00' 13" 

Dec 20 

2 X 1.5 

2 X 3.5 

Landolt 

98 

Ob''51™ 50" 

-00° 21' 17" 

Dec 21-24 

2 X 1.5 

2 X 3.5 


Table 1: CMOS observing log. The Landolt stars are photometric standards discussed in 
Landolt (1992). 
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bright foreground stars {r < 21), bright galaxies, and pixels which fall in the shadow of the 
on-instrument wavefront sensor. We also trim the outer perimeter of both helds where the 
SNR falls off due to the dithering process. The proportion of pixels lost as a function of 
galactocentric radius is shown in Fig. 6&, discussed below. The star counts are normalized 
to the dotted curves in this figure. 


5. Stellar photometry 

Fields 1 and 2 were analyzed using the IRAF daophot package (Stetson 1987; Davis 
1994). For each held, catalogues of Sa, da and 5a source detections were produced with the 
daoGnd routine. For each catalogue, preliminary aperture photometry was obtained using 
the phot task. The routine psfselect was used to select approximately the 100 brightest 
stars which were then inspected visually. Isolated stars (i.e. without bright neighbours) 
which were not saturated were chosen for psf stars. There were about 15 such stars in each 
held and they were used as input for the task psf which iteratively computes the point- 
spread function model for both helds. Finally, the task allstar was run to simultaneously ht 
this model to stars from the daohnd catalogues and determine their photometry. The task 
daophot identihed and returned htted parameters for 15671 and 4468 sources in Fields 1 and 
2 respectively; these numbers dropped to 15390 and 4291 sources after trimming the held 
edges. 

We examined the residual maps produced by allstar to evaluate psf subtraction with 
diherent thresholds. In general, the residuals were slightly cleaner for Field 1 when compared 
with Field 2 but for both helds the residuals were less than 1% of the total source hux. In 
Fig. 3, we show a small region from Field 1 before subtraction, after subtracting sources 
above 5a and after subtracting sources above 3a. The excellent photometry and source 
subtraction is evident. 

The photometric calibrations were established with two Landolt helds (see Table 1) 
for which there were half a dozen useful stars in each held (Landolt 1992). The equivalent 
magnitude for a single photoelectron was determined to be mp = 28.20, i.e. 0.1 mag fainter 
than what is quoted at the Gemini CMOS web site. The photometric standards showed that 
the photometry was better than 0.04 mag from night to night, and better than 0.02 mag in 
consecutive exposures. 

In Fig. 4, we show the magnitude distribution of detected sources found with daophot. 
Our catalogue is complete to r' = 27.0 (3a), r' = 26.7 (4a), r' = 26.3 (5a). For comparison. 



the expected CMOS performance can be assessed with the web-based calculator^, with fair 
agreement. Our calculation assumes a G5III spectrum. For 50 percentile observing condi¬ 
tions, an Rc = 26 point source in the r' hlter requires 8100 sec (9x 15 min dithered exposures) 
for a 5a detection. The conversion from to r' is given at the end of this section. 

Fig. 4 therefore establishes that the observations were undertaken in better than 50 
percentile conditions. As a further check on the integrity of our data, the star counts in the 
overlap region between Fields 1 and 2 were compared. We hnd 10% more stars in Field 1 
compared to Field 2, which is broadly consistent with the higher average background and 
slightly degraded seeing in the outer held. As we stressed earlier, the psf subtraction was 
cleaner for Field 1 compared to Field 2, and this may contribute to the difference. The 
magnitude offset between Field 1 and Field 2 was less than 0.04 mag and consistent with 
measurements from the Landolt standards. The uniformity between helds is evident from 
Fig. 4. 

It is convenient to establish a correction from r' to Rc for two reasons. Most surface 
brightness prohles published to date are measured in traditional Johnson-Cousins bands 
{e.g. Kregel et al. 2002). Secondly, in a later section, we correct the star counts for the 
background contribution due to distant galaxies. Again, galaxy count surveys to date are 
quoted in traditional bands {e.g. Metcalfe et al. 2001). Girardi et al. (2002; 2004) derive the 
requisite transformation from stellar atmospheric models. In the temperature range 4000- 
10,000 K, r' — Rc changes from about 0.3 to 0.15, with a stronger dependence for cool stars. 
This compares favourably with Fukugita et al. (1996) who hnd r' — Rc = 0.16(14 —-Rc) + 0.13 
such that for stars with 14 — Rc colours of 0.2 to 0.7, r' — Rc lies in the range 0.16 to 
0.24. Thus, to a good approximation, we can convert our magnitude system to Rc using 

r' - R^ = 0 . 2 . 


^See www.gemini.edu/sciops/instruments/instrumentlTCIndex.html 


Fig. 4.— r' magnitude distribution for all sources detected by daophot. The three curves 
are from the (top) 3a catalogue, (middle) 4a catalogue, (bottom) 5a catalogue. 
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6. Results 

Previous studies have shown that the slope of the radial light distribution of some disks 
becomes abruptly steeper at an outer break radius {e.g. Pohlen et al. 2002; Kregel et al. 
2003); the reason for this break is not known. We now seek to establish whether the stellar 
disk truncates rapidly beyond the break or extends faintly into the HI disk. Both outcomes 
have implications for how most of the baryons must have settled during dissipation. 

We derive the radial prohle of NGC 300 from the catalogue of htted stars. The centre 
of the galaxy was taken from the NASA/IPAC Extragalactic Database. Since the disk is 
inclined, we need to dehne elliptic annuli; each annulus has a radial (annular) thickness of 
100 pixels. The form of the annuli is illustrated in Fig. 2 for which we adopt a PA = 110° 
and inclination i = 42° (Kim et al. 2004). The angle subtended by our helds on the sky is 
33° with respect to the nucleus; this translates to 43° projected onto the plane of the disk, or 
about an eighth of the total azimuthal extent. Note that the disk major axis passes through 
the middle of Field 1. 

In Fig. 6a, we show the radial profile in NGC 300 derived from our source counts prior 
to removal of the background galaxy counts. The form of the prohle is essentially identical 
for the 3(7, 4a and 5a catalogues. The radial prohle is normalized by keeping track of the 
number of pixels in each annular bin. In Fig. 6b, we see that a typical bin has almost 200,000 
pixels. The top hat distributions are similar in form, and overlap at 18.5'. The combined 
distribution gives uniform coverage from about 12.5' out to 24.5'. The dotted lines indicate 
the small corrections required for pixels which are masked out due to bright foreground stars, 
bright galaxies and so on (see Fig. 2 and §4). The star counts are normalized to the dotted 
curves in this hgure. 

In Fig. 6a, there is a gentle transition at 19' (11.4 kpc) where the prohle appears to 
hatten oh. The transition is somewhat emphasized by the nature of the log-linear axes 
traditional for this held. This happens to fall in the overlap region between Fields 1 and 2, 
as we show in Fig. 6b. However, the change in slope is not an artefact of joining two data 
sets, since the form of the transition is identical in both helds. The transition is an artefact 
of a faint background galaxy population observed through the outer disk of NGC 300. 


Fig. 5.— The derived r' magnitude in daophot for all htted sources as a function of major 
axis radius. The decline in source counts as a function of radius is evident. Note the 
incompleteness at 12' due to crowding in the bright inner disk. Note also the excellent 
photometric consistency of the data between Fields 1 and 2. 
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In order to derive the intrinsic surface brightness prohle for NGC 300, we need to remove 
the background galaxy counts. In Fig. 7, we show three possible levels for background 
subtraction. In Fig. 7a, we take the average galaxy counts from the Hubble Deep Fields as 
derived by Metcalfe et al. (2001), where we have taken r' — Rc = 0.2 (see §5). Fig 7h takes 
the outer third of Field 2 to dehne the background counts. This is 10% higher than the 
averaged Hubble Deep Field counts. Given that the expected variance in a 5.5' held down 
to r' = 27 can be as high as 50% from held to held, this would not be unexpected. The 
diherence in counts between the Hubble Deep Fields when integrated down to r' = 27 is 
about 15% (N. Metcalfe, personal communication). 

For Figs. 7a and b, the key implication is that the disk extends to at least 24' (2. 2 /^ 25 ), 
close to 10 optical scale lengths for this low luminosity galaxy. It is possible to remove 
the transition at 18.5' completely with a background level which is 100% higher than the 
Hubble Deep Field counts. This outcome is presented in Fig. 7c. In contrast to Figs. 7a and 
b, the interpretation here would be very diherent. The disk now appears to truncate very 
sharply at 18', or about 50% further out in radius than the Holmberg radius. Such a cut-off 
corresponds to about 6 optical scale lengths, although note that the disk is still detected 
to 1.9i?25- While we consider this high background level to be unlikely, it cannot be ruled 
out until more of the outer perimeter of NGG 300 has been mapped under similar observing 
conditions. 

We are not able to separate faint stars from faint background galaxies using the GMOS 
images. We attempted to separate these with Sextractor but most of the sources identihed 
as galaxies could not be distinguished from stars convolved with the mildly elliptic psf. 
We therefore make the basic assumption that stars and galaxies have the same magnitude 
distribution in the range r' = 23 — 27. We can then directly convert the magnitudes and 
number counts in Fig. 7b to a mean stellar surface brightness as a function of radius. 

In Fig. 8, photometric data for NGG 300 are superimposed from three sources: (i) 
Garignan’s photographic measurements in Bj, (ii) I band measurements from Kim et al. 


Fig. 6.— (a) Source counts per arcmin^ as a function of radius in the outer GMOS helds 
for (top) 3a, (middle) 4a and (bottom) 5a detections. The radial binning is performed in 
projected elliptic annuli with 100 pix width (14.6"). The background galaxy counts have not 
been removed. The error bars are given by 1/\/N where N is the number of stars measured 
within the annular bin. (b) The upper solid histograms show the total number of GGD pixels 
per radial bin in Fields 1 and 2. The lower solid curves show the number of pixels rejected 
by the mask in Fig. 2. The dotted curves show the impact of the rejected pixels on two 
upper curves. 
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(2004), and (iii) our new CMOS r' measurements. The r' stellar surface brightness prohle is 
derived from the 3a catalogue presented in Fig 7b after background subtraction. Note that 
this prohle reaches exceedingly faint levels of = 30.5 mag arcsec“^ (at 2a when averaged 
over the last two bins in Fig 7b). 

We have offset Kim’s measurements downwards by 0.5 mag {r'—I = 0.5), and Carignan’s 
measurements upwards by 1.1 mag {Bj—r' = 1.1). Our population synthesis model in Fig. 10 
indicates that an appropriate correction is Bj — r' = 1.4 — 1.5. After applying our correction 
for stellar incompleteness below, we cannot account for the 0.3 mag offset with respect to 
Carignan’s data; this discrepancy was also noted by Kim et al. when they compared their 
data with Carignan (1985). 

In Fig. 9, the logarithmic axes emphasize the different htted components for the stellar 
core, the bulge and the disk. The htting procedure is described by Kim et al. (2004). 
Roughly speaking, the data conform to an exponential disk out to 24'. The substructure 
along the curve may be indicative of our restricted azimuthal coverage. The bulge cannot 
be strongly constrained in Kim’s htting. However, our hgure does make the point that there 
is no compelling evidence for a transition from disk to bulge or halo stars at the outermost 
observed extent. 

The fraction of missing light was estimated from a theoretical colour-magnitude diagram 
derived from the StarFISH code (Harris & Zaritsky 2001). In our models, we explore two 
possibilities: (i) an intermediate aged 8 Gyr stellar disk with [Fe/H]=-0.87, consistent with 
the metallicity derived from the red giant population (Tikhonov et al. 2005); (ii) an old 
12 Gyr disk with metallicity [Fe/H]=-0.7 comparable to the thick disk of the Galaxy. We 
specify a Salpeter initial mass function and generate 10^ stars for each model. The mass 
function uses a (minimum, maximum) stellar mass of (0.1,100 ) Mq, and a binary fraction of 
25%. The adopted isochrones are taken from http://www.te.astro.it/BASTI (Pietrinferni et 
al. 2004). 

The resulting colour-magnitude diagrams are shown in Fig. 10. Note that the models 
predict that the CMOS data recover (i) 54.6% (ii) 55.8% of the total light at our 3a threshold. 
Therefore, we estimate our completeness at 55% in the r' images. This would indicate that 
our surface brightness prohle in Figs. 8 and 9 should be shifted upwards by -0.65 mag in 
order to correct for the missing light. (We do not attempt to correct for the internal dust 
extinction.) However, we note that the galaxy inclination [i = 42°) renders the surface 
brightness prohle brighter by about -0.32 mag. Therefore, the true upward correction is 
expected to be -0.33 mag. 

With our derived conversion from r' to Rc (§6), this indicates that the CMOS data reach 
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down to an effective surface brightness of = 30.3 mag arcsec”^. The CMD simulations 
in Fig. 10 predict an outer disk colour B — I ^ 1.7. The equivalent surface brightness in B 
reached by the CMOS data is at the level of /i^ = 32 mag arcsec“^. These are levels which 
cannot be reached reliably in diffuse light imaging as we discuss in the Appendix. 


7. Discussion 

Our data clearly reveal the presence of a faint stellar disk out to at least a radius of 24' 
(14.4 kpc) or about 10 disk scale lengths. Our data extend far beyond the surface photometry 
of Kim et al. (2004) and show that the disk remains exponential to the radial limits of our 
counts. There is no evidence for a break in our data, which cover an angular wedge of 
about 45° within the deprojected disk - the stellar surface brightness is well described by 
an exponential over almost 10 scale lengths. There is no evidence of a bulge or a halo 
component. 

Indeed, in three wide-held photometric studies of NGC 300, there is no compelling 
evidence for a ‘break radius’ in the galaxy luminosity prohle (de Vaucouleurs & Page 1962; 
Carignan 1985; Kim et al. 2004). Carignan (1985; Fig. 5) traces the light in a deep UK 
Schmidt plate out to [xb = 27.5 mag arcsec”^ at r = 12.5'. While there may be a hint at 
a break in Carignan’s Fig. 14 at about 10', his earlier hgures demonstrate quite different 
structure when the galaxy is split into two halves. 

Since our stellar counts extend to the outer HI disk, it is important to consider the 
possible effects of a warped stellar disk. Puche et al. (1990) hud evidence for a warp in the 

Fig. 7.— Star counts as a function of major axis radius for the 3a detections. The counts 
in the overlap regions have been averaged. We show the counts after subtracting 3 different 
background galaxy count levels indicated by the horizontal line: (a) predicted counts from 
averaged Hubble Deep Fields; (b) average source counts in the outer 8 radial bins (annular 
width of 800 pixels) assuming that all sources are galaxies; (c) an artihcially high background 
from assuming that all detections in the range 20'-24' are background galaxies. In each case, 
the background level is shown by a solid line. The error bars are discussed in the caption 
for Fig. 6. 
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Fig. 8.— The observed r' surface brightness prohle of NGC 300 derived directly from the 
star counts. The inner crosses are taken from Kim et al. (2004): these are I band measure¬ 
ments shifted downwards by 0.5 mag. The open circles are taken from the Bj photograpic 
measurements of Carignan (1985) and are shifted upwards by 1.1 mag. The CMOS data 
points are derived from the 3a source catalogue presented in Fig. 7b. After correcting for 
incompleteness (—0.65 mag) and inclination (-1-0.33 mag), we obtain the intrinsic surface 
brightness prohle (see §6); this amounts to shifting all points upwards by 0.32 mag. 


Fig. 9.— The r' surface brightness prohle of NGC 300 presented in log radius to emphasize 
the contribution of the diherent htted components. The data points are described in Fig. 8. 
The continuous curves are the diherent htted components, for which the functional forms 
are given in Kim et al. (2004). The model for the core is not shown. 


Fig. 10.— The colour-magnitude diagram for a stellar population with (left) mean age 8 Gyr 
and [Fe/H]=-0.9, (right) mean age 12 Gyr and [Fe/H]=-0.7. The intrinsic stellar brightness 
reached by our data is indicated by the horizontal dashed lines for la, 3a and 5a detections. 
The right hand hgure of each panel shows the fraction of received light from the stellar 
population above our threshold level. The greyscale indicates the number of stars used in 
the simulation. 
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HI disk beyond 10' but this does not necessarily indicate that the stellar disk is warped in 
the same way, if at all (e.g. Garcia-Ruiz, Sancisi & Kuijken 2002). Dramatic stellar warps, 
which are rarely observed, are most often associated with strong galaxy interactions. There 
is no obvious interacting companion observed in the vicinity of NGC 300. A weak stellar 
warp may exist and its effect would be to slow the radial decline in the stellar luminosity 
prohle. However, there is no such effect evident in either Fig. 7 or Fig. 9. 

The surface brightness of the outer disk of NGG 300 extends smoothly down to at least 
30 r' mag arcsec”^ which corresponds to a stellar surface density of order 0.01 Mq pc“^. To 
the best of our knowledge, such an extended and very diffuse population of stars has not 
been seen before. What is its origin? 

A key issue is just how far the HI disk extends. Puche et al. (1990) have mapped the 
outer disk in NGG 300 and detect HI out to 32', beyond the extent of the newly established 
optical disk. The HI surface density prohle is roughly constant at 0.5 — 1 Mq pc“^ in the 
outer disk. The surface density prohle of the optical disk falls with radius, and is roughly 
comparable to the HI at the Holmberg radius, assuming M/= 1. Beyond here, the optical 
surface density falls by a factor of 100 at a radius of 20' or more. 

The existence of stars at such low surface densities is difficult to understand in the 
conventional picture. How does a spiral density wave initiate and propagate - if this is 
indeed what triggered star formation in the outer disk - in such a rarehed medium? 

The Q criterion is often used to establish whether a disk is unstable to axisymmetric 
modes, viz. 

^ ~ 3.36GS 

where a is the internal dispersion of the HI gas which stabilizes the disk on small scales 
(up to the radius of the Jean’s mass). On larger scales, the disk is stabilized by differential 
rotation which is embodied by k, the epicyclic frequency of the disk. In contrast, the disk 
is destabilized by its own surface density, S. The disk is locally stable against axisymmetric 
modes if Q is substantially above unity. 

We derive /t(r) and S(r) from the data of Puche et al. (1990) and assume <7 = 5 km 
s“^ in the outer disk. We thus hnd that Q is everywhere greater than 5±2 beyond 10' which 
would appear to argue against star formation and the existence of spiral density waves here. 

However, the observations of Puche et al. (1990) not only show that the HI extends 
well beyond the extent of our stellar disk but also reveal well-dehned, tightly wound spiral 
structure in the HI which is particularly prominent in the SE at 20'. While this spiral 
structure is difficult to explain, it is interesting that HI spiral arms are increasingly observed 
in the outer parts of disk galaxies (TA Oosterloo, personal communication; Bureau et al. 
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1999; Quillen & Pickering 1997). Bureau et al. found a similar problem for the HI disk of 
NGC 2915, which also shows marked spiral structure at very large radius; they invoked the 
effect of a tumbling triaxial halo potential to drive the spiral structure in this apparently 
stable outer disk. 

It is known from photometric imaging studies (e.g. de Jong 1996) that most Hubble 
types have colour gradients from the centre outwards which arise from (i) the decreasing dust 
content; (ii) declining stellar abundances; and (hi) changing stellar populations due to radial 
variations in star formation history. The general trend is that disks become progressively 
bluer at larger radii and this has been ascribed to the stellar population becoming younger 
on average and progressively more metal poor (de Jong 1996). 

However, the inference that the stellar populations become progressively younger with 
radius contrasts with what we know from recent stellar photometry of the outermost reaches 
of nearby galaxies. These studies incorporate high quality stellar photometry and show, by 
comparison with well calibrated globular cluster data, that the outer giant populations are 
metal poor and old. Recent examples include the outer disks of M33 (Galleti, Bellazzini & 
Ferraro 2004) and M31 (Ferguson & Johnson 2001). 

The data for NGG 300 are limited at present but Tikhonov et al. (2005) hnd that metal 
poor ([Fe/H]=-0.87) giants dominate the light at large radius and these are most likely to be 
old. We note that Tikhonov et al. attribute various regions of NGG 300 to the thick disk or 
to the halo which is difficult to do in a relatively face-on galaxy. One can really only rely on 
indirect arguments based on our Galaxy, but these are very unreliable. For example, their 
outermost held (S3) corresponds to a deprojected radius of about 11 kpc or roughly 6 scale 
lengths. In our Galaxy, this corresponds to about 20 kpc where we know nothing about the 
stellar disk at the present time. 

So is the outer disk in NGG 300 young or old? When one considers the effects of changing 
extinction and metallicity, a moderately blue outer disk does not necessarily indicate a 
younger stellar population. (Galleti et al. (2004) discuss how age and metallicity effects can 
be disentangled in outer disk/halo observations.) Within the context of cold dark matter 
models, the answer to this question has key implications for the formation scenario for the 
outer disk (e.g. Ferguson & Johnson 2001). 

If follow up measurements reveal that the outer disk is moderately young, we would 
expect to hnd evidence of spiral density waves in the stellar population even at the implied 
large radii (>12 kpc), although these spiral arms may be moderately dihuse and hard to 
recognize. As mentioned above, there is clear evidence of spirality in the outer HI disk. If 
stars are forming today, it may be that the outer HI disk largely comprises cold, compact 
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clouds with internal dispersions lower ( ^ 1 km s“^) than one normally associates with 
beam-averaged observations of outer disks (~ 5 km s“^). If the HI gas was to cool to 100 
K, and nonthermal (e.g. turbulent) velocities were ^ 1 km s“^, the Q criterion could be as 
low as unity. Such clouds would be gravitationally unstable today and the outer disk would 
have a moderately young stellar component. But we stress that, to our knowledge, the HI 
velocity dispersion has not been measured in distinct clouds in the outer disk of NGC 300. 

If new measurements reveal that the outer disk is old, then the current high value of Q 
here may simply reflect that the HI has been largely used up today, but that there was a time 
~10 Gyr ago {z ~ 1) when there was far more gas than we observe today. If our assumption 
of the cloud dispersion velocity still holds, this would have the effect of raising the surface 
density of the disk, and lowering the Q parameter to the point where local instabilities induce 
star formation. 

Hydrogen in the outermost reaches of galaxies may be evidence for gas accretion, either 
from gas settling onto the dark halo (cold accretion), or from a galaxy merger event (hot 
accretion). If the material in the slow warp in NGG 300 is the result of gas accretion to the 
outer disk, it is plausible that the settling process produced low levels of star formation. An 
example of such a process is dramatically illustrated by the gas disk settling within the giant 
elliptical NGG 5128 where a young bine popnlation is clearly evident (Graham 1979; Bland, 
Taylor & Atherton 1987). 

But how does an accretion process sncceed in maintaining the continnity of the exponen¬ 
tial disk? This is particularly puzzling when the surface brightness distribntion is contrasted 
with that of most other spirals, which show a clear break in their light distribution at 3 to 5 
radial scale lengths. Gould the Lin-Pringle viscous evolution of a star-forming disk explain 
the persistence of the exponential disk of NGG 300 to ten or more scale lengths? 

Another possibility is that the stars were scattered from the inner disk. N-body simula¬ 
tions of disk evolution indicate that radial mixing is strong (Sellwood 2002). A single spiral 
wave near co-rotation can pertnrb the angular momentum of a star by ~20% moving it either 
inwards or outwards by several kiloparsecs (Sellwood & Kosowsky 2002). This process could 
help to wash out any stellar disk edge bnt it is difficult to see how this explains the very 
extensive stellar disk, with a simple exponential density prohle, observed here. 


JBH would like to dedicate this work to his brother Simon who died in December 2004. 
We thank the staff of the Gemini Observatory for their tireless efforts and their pursnit of 
excellence. MV was snpported by an AAO Stndent Fellowship in the austral winter of 2004. 
BTD was snpported in part by NSF grant AST-9988126. We are grateful to Laura Dunn 
for her help with running the StarFISH code, and to the referee for a carefnl examination of 
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A. Starlight scattered by dust in the outer disk 


In our study, we reach an effective surface brightness of = 30.5 mag arcsec”^ (2o'). This 
level is exceedingly faint in the sense that the equivalent surface brightness in i? or is fainter 
than 32 mag arcsec“^ for an old stellar population. This raises an important question: Is 
there a practical limit to galaxy photometry? Here, we demonstrate that in principle the 
metal-poor outer HI envelope has a sufficient column of dust to scatter detectable levels of 
starlight from the optical disk. Thus, we propose that only star counts can reliably trace 
luminosity profiles in galaxies at this faint level or even fainter. 


Let us consider a disk with inclination angle i. Starlight from the inner regions of the galaxy 
will illuminate the dust grains, resulting in scattering. Consider a region in the galaxy where 
the H nucleon column density normal to the disk is Ah. If we approximate the illuminating 
starlight as due to a point source at a distance D with luminosity located at the galactic 
center, then the dust grains at a radial distance r from the galactic center will result in 
scattered intensity ly, where 


^ ^ Ly Ah 
47 rr^ cos i 




The scattering phase function can be written 


(Al) 


-^x(A, &) = y^ — 
^ uh 



(A2) 


for which {d(7/dQ)j is the differential scattering cross section for grain type j, at wavelength 
A and scattering angle 0*, and Uj/n-n is the number of grains of type j per H nucleon. Eq. 
(Al) can be rewritten 

I- ^ J_ \ ( Fx(A.e.) \ (A3) 

29 mag arcsec ^ cosf yarcsec/ \10^®cm ylO ^"^cm^sr ^y 

where mgai is the apparent magnitude of the galaxy. For NGC 300, the NED data base gives 
^^gai = 8.95 in the B band rising to mgai = 7.5 in the R band. 


Fig. 11.— F(A, 0) or Milky Way dust with Ry = 3.1 (lower panel) and SMC bar dust with 
Rv = 2.9 (lower panel) at selected wavelengths A (SDSS z’, i’, r’, g’, and u’, and Cousins ff, 
Rc, and Vc), as a function of the scattering angle 0*. 


Fig. 12.— Degree of polarization for scattering by Milky Way dust with Ry = 3.1, and SMC 
bar dust with Ry = 2.9 at selected wavelengths A, as a function of the scattering angle 0*. 
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If the location in the galaxy is specihed by the observed position angle ip (measured from the 
minor axis, with = 0 corresponding to the side of the galaxy nearer to us) and displacement 
0 from the center, then 


_ h __ cosi(l + tan^^) / 0 \ ^ / iVn ^ / Fx{\Qs) \ 

29 mag arcsec”^ 1 + cos^ i tan^ yarcsec/ ylO^^cm^^y y 10“^^ cm^ sr“^/ 

The scattering angle 0^ is related to inclination i and position angle tp by 


XQ0-4(14-mgai) 


(A4) 


leading to 0^ 


Qs = arccos 


sini 

y/l + cos^ i tan^ tp 


90° for a face-on galaxy (i = 0). 


(AS) 


Weingartner & Draine (2001) have developed dust models that reproduce the observed 
wavelength-dependent extinction in the Milky Way, LMC, and SMC. These models are also 
consistent with the observed infrared and far-infrared emission (Li & Draine 2001, 2002). 
The scattering properties of these dust models have been calculated by Draine (2003). Fig. 
11 shows F{Qs) from Draine (2003) at selected wavelengths for MW dust and SMC dust. 
Note that F{\, Qg) for the SMC is lower than for the Milky Way by approximately a factor 
10, which is approximately the metallicity of the SMC relative to solar. Thus, it seems 
reasonable to estimate Fx (A, 0s) ~ (Z/Zq)Fmw(A, 0s) as an estimate for Fx in another 
galaxy. 

Puche et al. (1990) measured iV(HI) ^ 9 x 10^® cm“^ in NGC 300 at 1000". For the outer 
disk of NGC 300, the red giant colour distribution indicates [Fe/H] = -0.87 (Tikhonov, 
Galazutdinova & Drozdovsky 2005); we take ZjZQ ^ 0.2 for the gas. Evaluating eq. (A3) 
for a point on the major axis 1000" from the center, we hnd / ~ 31.5 mag arcsec“^, or 
about 4% of the stellar surface brightness at this point. Scattered starlight makes a small 
contribution to the total surface brightness for NGC 300 at r = 1000", but for other cases 
(e.g., galaxies where the stellar surface density declines more rapidly) scattered starlight 
could contribute a larger fraction of the total surface brightness. 


Scattered starlight, if present, will be substantially polarized. In Fig. 12, we show the degree 
of linear polarization expected for the scattered light as a function of scattering angle, where 
we assume that the galactic starlight incident on the dust grains is unpolarized. For a 90° 
scattering angle (e.g. for a face-on galaxy) the polarization is expected to be ~ 45% in red 
light, with smaller values at shorter wavelengths. 


In addition to scattering light, dust grains appear to luminesce in the red when illuminated 
by ultraviolet light (see Witt & Vijh 2004, and references therein). This so-called “Extended 
Red Emission” can be comparable in intensity to the scattered starlight in the red. 
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